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Abstract 

We solve the monomer-dimer problem on a non-bipartite lattice, the simple quartic lattice with 
cylindrical boundary conditions, with a single monomer residing on the boundary. Due to the 
non-bipartite nature of the lattice, the well-known method of a Temperley bijection of solving 
single-monomer problems cannot be used. In this paper we derive the solution by mapping the 
problem onto one on close-packed dimers on a related lattice. Finite-size analysis of the solution is 
carried out. We find from asymptotic expansions of the free energy that the central charge in the 
logarithmic conformal field theory assumes the value c = — 2. 

PACS numbers: 05.50. +q,02. 10. Ox,11.25.Hf 
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I. INTRODUCTION 



An outstanding unsolved problem in lattice statistics is the monomer-dimer problem. In 
this problem diatomic molecules adsorbed on a surface are modeled as rigid dimers occupying 
two adjacent sites and lattice sites not covered by dimers are regarded as occupied by 
monomers. While the case of pure dimers has been solved in 1961 by Kasteleyn [1] and 
by Fisher and Temperley [2, 3], the general monomer-dimer problem has proven to be 
computationally intractable [4]. 

In 1974, Temperley [5] introduced an intriguing bijection mapping the dimer problem 
with a single monomer at the corner of a finite M x N lattice to the counting problem 
of spanning trees on a related lattice, thereby providing an alternate way of deducing the 
solution. The method of Temperley bijection has since been extended to the case when the 
monomer resides on other specific boundary sites [6]. However, the success of the Temperley 
bijection apparently relies on the fact that the lattices being bipartite; it does not work for 
non-bipartite lattices. In this paper, we consider one non-bipartite lattice, a rectangular 
lattice with a cylindrical boundary condition. By using an alternate mapping formulated 
recently by one of us [7, 8], we solve the monomer-dimer problem on this lattice when a 
single monomer resides on the boundary. We also clarify the mathematical content of the 
solution by carrying out finite-size analysis of the solution. 

II. SINGLE MONOMER ON THE BOUNDARY OF A CYLINDER 

Consider a simple quartic lattice C consisting of an array of N rows and M columns 
embedded on the surface of a cylinder with periodic boundary conditions imposed in the 
horizontal direction. See Fig. 1(a) for an illustration. For MN odd, hence both M, N 
odd, the lattice is not bipartite. But the lattice can be fully covered by one monomer and 
(MN — l)/2 dimers. We consider the problem of evaluating its generating function when 
the single monomer resides on the boundary. 

On first sight, one would attempt to use the Temperley bijection of mapping. However, 
it can be readily verified that the attempt invariably fails, apparently due to the fact that 
£ is not bipartite. Instead, we adopt an alternate formulation devised by one of us [7, 8] 
which does not make use of the Temperley bijection. 
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FIG. 1: (a) A simple quartic lattice C consisting of an array of N = 3 rows and M = 5 columns 
embedded on the surface of a cylinder, (b) A self-dual lattice C derived from C by adding a new 
site S connected to all M sites of one boundary of C. (c) A oriented lattice C" constructed from 
£ by keeping only one edge connecting to S. A phase factor i is associated to all x dimers. 

Denote the desired generating function by 

G MD (x,y)= J2 xni y n2 > (!) 

config 

where the summation runs over all monomer-dimer configurations with a single monomer 
on one of the two boundaries, x > and y > are the weights of, respectively, horizontal 
and vertical dimers as indicated in Fig. 1(a), and ri\ and 112 the numbers of horizontal and 
vertical dimers subject to n 1 + n 2 = (MN — l)/2. For quick reference we first give the final 
result which holds for M, N > 3, 

M-1 JV-1 

G MD (x, y) = 2Mx^-^ N -^ f[ II ( Ax2 sin2 2 ~jf + V co§2 ATTl) • ( 2 ) 

m=l ri=l ^ ' 

In contrast, the monomer-dimer generating function with a single monomer on the boundary 
of an M x N net with free (open) boundaries is [6] 

M-1 Af-1 

G f ^(x, y) — (M + N — 2 )*<"-i)/y"-i>/* f[ f[ Lx> cos 2 ^ + V cos 2 , 

m=l n=l ^ ' 

(3) 

where the factor M + N — 2 is the number of equivalent boundary sites where the monomer 
can reside. 
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Results of enumerations of (2) and (3) for small lattices are shown in Table I. 



TABLE I: Enumerations of monomer-dimer configurations. 



M x N lattice 


Gmd(1, 1) given by (2) 


G md(M) given by (3) 


5x5 


3,190 


1,536 


5x7 


53,010 


24,150 


7x5 


56,434 


24,150 


7x7 


3,118,178 


1,204,224 


7x9 


171,527,426 


57,961,134 


9x7 


165,771,810 


57,961,134 


9x9 


29,845,632,402 


8,921,088,000 



To derive (2) we consider first the close-packed dimer problem on a related lattice £ 
constructed from C by connecting all M sites on one boundary to a single new site S as 
shown in Fig. 1(b). Dimers connecting boundary sites to S all carry weight 1. It is of interest 
to note that the lattice £ is self-dual and that the lattice has been considered previously by 
Lu and Wu [9] in the context of Ising partition function zeroes. 

Denote the generating function of close-packed dimers on £ by G £,{£', x, y). Since in a 
close-packed configuration S must be covered by a dimer (of weight 1), and the dimer must 
end at one of the M equivalent boundary sites which can be regarded as being occupied 
by a monomer on C, there exists a correspondence between dimer configurations on £ and 
monomer-dimer configurations on £ We are led to the identity 

G M D(x,y) = 2G D (£;x,y), (4) 

where the extra factor 2 comes from the fact that there are 2 boundaries on a cylinder. 

To evaluate G D (£;x,y) we introduce the lattice £ shown in Fig. 1(c) where S is con- 
nected to only one boundary site. Denote the generating function of close-packed dimers on 
£' by G D (C"; x, y). It is clear that we have the further identity 

G D (£;x,y) = MG D (£;x,y). (5) 

It remains to evaluate G D (£';x,y). But this is the problem solved in [7, 8]. 
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In the analysis given in [7], close-packed dimers on a lattice similar to C" are enumerated 
using the Kasteleyn approach [1]. Since our procedure follows closely that discussed in [7], 
we give an outline and highlight the difference. 

Orient edges of C" and associate a phase factor % to all x edges as shown in Fig. 1(c). The 
only thing new from [7] is that we need to ascertain signs of all terms in the Pfaffian are the 
same. However, it can be shown [10, 11] that this always is the case for M = odd. Then the 
desired generating function Gd(C";x,i/) is given in terms of the Pfaffian of a matrix A' [8], 



iS M -^ 2 G D (C";x,y) = Pf(A') = VfetA. 



(6) 



Here A' is the antisymmetric Kasteleyn matrix of dimension (MN + 1) x (MN + 1) for the 
lattice C" explicitly given by 

/ ••• 1 ••• \ 





A' 





-1 




(7) 



V o / 

where A is the Kasteleyn matrix of dimension MN x MN for C The position of the 
elements ±1 in the first row and column is that of the site {m, 1} connected to S (see 
below). Explicitly, A is given by 



A = ixS M <8> In + ylu ® T N , 
with 1m is the M x M identity matrix, Sm is the periodic M x M matrix 

/ 1 ••• -l\ 
-1 1 ••• 
-1 ••• 



(8) 



0-- 
\ 1 •• 



1 
-1 J 



(9) 



and T/v is the N x N matrix 



/ 1 
-1 1 
0-10 



■ N 



■ o o\ 

• 

• 



(10) 



••• 1 
\ ••• -1 0/ 

Note that we have Tm instead of Sm in the corresponding expression in [7]. 

Label elements of A by {to, n; to', n'}, where (to, n) specifies the column and row of the 
position a site. The determinant of the Kasteleyn matrix A' can be computed by Laplace 
expanding along the first row and first column leading to 



detA' = C(A; {to,1;to,1}), 



(11) 



where C(A; {to, 1; to, 1}) is the cofactor of the {to, 1; to, 1} element of A, and we have spec- 
ified the site connecting to S in Fig. 1(c) as {to, 1}. 

Since the cofactor C(A; {to, 1; m, 1}) is proportional to the product of the nonzero eigen- 
values of the matrix A, we need to determine the eigenvalues of A. This is done in the next 
section. 

III. EIGENVALUES OF THE KASTELEYN MATRIX A 



The matrix Sm can be diagonalized by the similarity transformation 

where Vm and its inverse V M X are M x M unitary matrices with elements 

1 



V M (m 1 ,m 2 ) 



i2mim,2ir/M 



M 



V M \ mi ,m 2 ) = 1 g-aminWM x < {mijm2 } < M , 
v M 

and VLm is an M x M diagonal matrix with the eigenvalues a; m of Sm as entries, 

2TO7T 

ui m = 2i sm — — , 1 < to < M . 

M ~ ~ 



(12) 



Similarly, as in [7], the matrix T/v is diagonalized by the similarity transformation 

U^TnUn = Tat, 

where Un and its inverse U^ 1 are N x N unitary matrices with elements 

tt l \ I ^ . ni . ( n x n 2 -n\ 

for 1 < {ni,n 2 } < iV, and r^r is an iV x N diagonal matrix having eigenvalues 7„ of T N as 
entries, 

nix 

7„ = 22 cos , 1 < n < N . 

Thus, the MN x MN matrix A can be diagonalized by the similarity transformation gen- 
erated by Umn — Vm <S>Un, leading to 

Um N AUmn = Amtv, (14) 

where A M n is an MN x MN diagonal matrix having eigenvalues \ m . a of A as entries, 

ix sin + y cos ^. J, 1 < m < M, l<n<N. 

Note that \ m , n vanishes at m = M, n = (N + l)/2. Elements of Umn and its inverse U M l N 
are 

U M N{m 1 ,n 1 ;m 2 ,n 2 } = V M (mi, m 2 )C/jv("i, ^2) 

^Miv{ m i> n i; m 2,W2} = V^ 1 (mi,m 2 )C/^ 1 (ni,n 2 ). (15) 
Using the identities sin(27r — 0) = — sin 9 and cos(7r — 0) = — cos 9 , the product 

M N-l 

p=n n v«, (is) 

m=l n=0 

(m,n)#(M,^i) 

where the product excludes the zero eigenvalue at (m,n) = (M, ^f-), can be rearranged as 

M-l JV-1 
m,=l n=l ^ 
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where the factor Q collects all factors with either n — (N + l)/2 or m = M, namely, 



M-1 

9 9 2to7t\ 
- sin 



^ = n - 4x 



m=l 



M 



II % 

/ n=l V 



• 2 cos 2 



nix 



N + 1 



( 1)( m-i)/2 ^(^+i: 



(18) 



after using the identities 

2mn 



N-l 
2 



m=l V 7 n=l V 7 

The expressions (17) and (18) apply to M, N > 3 and will be used in the next section. 



IV. EVALUATION OF THE GENERATING FUNCTION (1) 

We now compute the generating function (1). 

Combining (4)- (6) with (11), we obtain the following expression, 

Gmd (x,y) = 2M ^- M )/ 2 y/C(A;{m,l;m,l}) . 



(19) 



where C(A; {m, 1; to, 1}) is the cofactor of the (to, 1; to, 1) element of the matrix A. 

The computation of cofactors of a singular matrix like A requires special attention since 
the matrix does not possess an inverse. The difficulty was resolved in [7] by perturbing 
the matrix A slightly rendering it non-singular to permit an inverse. By carrying out this 
analysis details of which we refer to [7], one finds the cofactor 



C(A; {to, n; to', n'}) = 



U M n (rri, n; M, J~ 1 ) U M X N (M, N ^~ 1 ; to, n) 



P, (20) 



where U M n is the matrix diagonalizing A. Note that the index {M, ^-^-} is that of the zero 
eigenvalue. 

Elements of Umn and U^ N are given in (15). After combining with (12) and (13), we 
obtain from (20) 

2i n '- n 



C(A; {to, n; to', n'}) 
valid for general to, n, to', n'. 



rnr nir 

MiNTT) 8111 T sm — 



(21) 
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Finally, we combine (4)-(6) with (11) and (21) at {m! = m,n' = n = 1}, and arrive at 
the expression 

G MD(x ,y) = 2Mfi-^ ] j J ^ Ty (22) 

This yields the generating function (2) given in Sec. II after substituting with P given by 
(17) and Q by (18). We note that the result is independent of m as it should. 
Then, with the help of the relations 

N-l N-l M-1 M-1 

n=l v 7 n=l v 7 m=l v 7 m=l 

valid for any function F(-), the generating function (2) can be written in the equivalent 
form, 



M-1 Af-1 
2 2 



G MD (x, y) = 2Mx( M - 1 »V N - 1)/2 f[ fl ( Ax2 sin2 ^ + V sin2 ]v^Tl) ' M ' ^ = ° dd - 

m=l n=l ^ 7 

(23) 

It is convenient at this point to introduce a function 

■M-1 JV-1 -i 1/2 



H(z; M, N) 



n n 



m=0 n=0 

(m,n)^(0,0) 



= z MN - 1 H(l/z;N : M), any M,JV>1. (24) 
It will be shown in Appendix A that we have 

G M d(x, y) = R M , N (y, z) ^H(z; M, N + I), M,N = odd, (25) 
where z = x/y and 

MvIy MN ~ x 
^^'^ 2 = (JV + l)*"^*)' 

Sm(^) = sinh(Msinh- 1 (l/z)). 

The advantage of using (25) instead of (23) for the generating function is that the factor 
Rm,n{u,z) sorts out major contributions in the asymptotic expansions of the free energy 
(30) and (31) discussed below. 

Two equivalent expressions of H(z; M, N + 1) can be obtained by taking one of the 
products in (24) in a closed form. Taking the product over n, we obtain 

M-1 

H(z; M, N + 1) = (N + 1) J] 2 sinh [(N + l)u z (ff)] (26) 

m=l 



where 

u z (k) = sink' 1 (z sink) (27) 

is the lattice dispersion relation, and we have used the identities (A2) and (A4). 

Similarly, taking the product over m and making use of (A4) and the equivalence (24), 
we obtain 

N 

H{z; M,N + 1) — Mz M ( N+1) ~ 1 J J 2 sinh [Mu 1/z (^)] ■ (28) 

V. FINITE-SIZE ANALYSIS AND ASYMPTOTIC EXPANSIONS 

Define the "free energy" of the monomer-dimer system as 

F MtN (x,y) = -]nG M D(x,y) 

= -\nR M , N (y,z)-^lnH(z;M : N + l), (29) 

where we have made use of (24). We note that other than an overall factor [4 sin 2 (cot/M) + 
4 sin 2 (/3tc/N)~\, the function H(z;M,N + 1) is the special case of a = (5 = of a more 
generally defined function Z a ${z\ M, N + 1) introduced, and analyzed in details in [12, 13]. 
This permits us to use results of [12, 13] to write down a general expression for F M;N (x,y), 
which we shall not reproduce. Instead, we focus on the free energies 

F M = lim —F MN (x,y) and F N = lim ^-F MN (x,y) 

of infinite "strips" and their asymptotic expansions. 

The asymptotic expansions can be deduced by applying the Euler-MacLaurin summation 
identity to In H(z; M, N + 1). Using H(z; M.N + 1) given by (26) and (28), respectively, we 
obtain from (29) using (26) and (28), respectively, 

M-l 

2 ' ' 2 



M 1 M_1 



m=l 

oo 



M f | V /7T\2P-1 d 2p _ 2 ( Z ) (B 2p \ 

- M/ bulk + ^ ^ - j ^— — ^— J 



HZ ( ' 

= Mf hnlk + — I — )+... , (infinite length) , (30) 



10 



N 1 1 N 

F N = --ln(yz) + -smh- 1 (l/z)--J2"i/z{jfi) 



Nfhu\k + 2 /surface + 



P =i 



' TT \ 2P 1 rf 2 p_2(l/^) / B 2p 



N+lJ (2p-2)! V 2p 



= N/ bulk + 2/ surfacc + ^- f — ^ - ) (infinite poriincntor). (31) 



where 



1 1 r 

/bulk = --lny-— / uj z (k)dk 

2 2ir J 

1 1 Z" 7 " 

= --ln(^)-— / u 1/z (k)dk, 

1 If 71 
/surface = - sinh -1 (1/^) - — y u 1/z (k)dk, 

d 2p (z) are the coefficients in the Taylor expansion 



with d (z) = z, d 2 (z) = ~z(l+z 2 )/3, d 4 (z) = z(l + z 2 )(l + 9z 2 )/5, ■ ■ ■ , and B 2 = 1/6, B 4 = 
— 1/30, B 6 = 1/42, • • • , are the Bernoulli numbers. The two equivalent expressions of / bulk 
are obtained from (30) and (31), respectively. We remark that the equivalence of the two 
expressions is verified by the intriguing integral identity 

-f [ siniT^zsinfl) - smli" 1 (± sin #)] dQ = \nz (32) 
71 Jo 

obtained by noting that the derivative of the left-hand side of (32) with respect to z reduces 
to \ j z after carrying out the integration. 

The general theory of finite-size analysis [14-16] dictates that the free energy per unit 
length of a lattice model at criticality on an infinitely long strip of width M assumes the 
form [16] 

Fjsf = M /bulk + /surface + + ' " , (33) 

in an asymptotic expansion where /bulk and /surface are free energy densities of the order of 
0(1) and A is a constant. Unlike the free energy densities, the constant A is universal and 
its value is related to the central charge c in the logarithmic conformal field theory in a 
relation which depends on the boundary conditions in the transversal direction. Explicitly, 
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A is proportional to an effective central charge c c g = c — 24 h m i n , where c is the central 
charge characterizing the universality class of the lattice model, as [14, 17] 

A = — — c c ff = — — (c — 24 h min ) on a cylinder of infinite length, (34) 
A = — — c c fj = — ^( c ~~ 24/i min ) on a cylinder of infinite perimeter, (35) 

where the number h m i n is the smallest conformal weight in the spectrum of the Hamiltonian 
with the given boundary conditions and ( is an anisotropy factor. In our case we find from 
(30) and (31) that ( — z and \/z, and A = it z/ 12 and tc/12z, respectively, in (34) and (35). 

To retain the characteristics of a monomer on the surface, we consider a cylinder of infinite 
perimeter in a geometry which retains two surfaces. Therefore we use (35) and (31), or Fn, 
for which the boundary condition in the transverse direction is free (open) boundaries. It is 
known [18] that for free (open) boundaries h m i n = 0. Hence we deduce the central charges 

c = c cff = -2. (36) 

On the other hand, if one uses (30), or Fm, the system is an infinitely long cylinder with 
a perimeter M. The two physical boundaries of the lattice are located at infinity so the 
existence of a monomer on the boundary is immaterial. The situation reduces to that of a 
pure dimer problem studied in [17]. For M = odd we are considering, the analysis of [17] 
also gives A = n(/12 as in (31). However, for M odd, the boundary in the transverse 
direction is "frustrated" requiring special attention. It is argued in [17] that in this case one 
should use (35) with /i m ; n = 0. This again leads to the same central charges (36). 

We remark that the c = — 2 central charge has been reported previously [6] in the solution 
(3) of a single monomer on the surface of a rectangular net with free (open) boundaries. 

VI. SUMMARY 

We have derived the closed-form expression of the monomer-dimer generating function 
for a non-bipartite rectangular lattice under cylindrical boundary conditions with a single 
monomer confined to reside on the boundary. We have also carried out a finite-size analysis 
of the free energy. Asymptotic expansions of the free energy of strips of infinite lengths in 
the periodic and free (open) directions are obtained using the Euler-MacLaurin summation 
formula. We find the central charge in the framework of the logarithmic conformal field 
theory to be c = — 2. 
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Appendix A 



In this appendix we establish the expression (25) for the generating function. 
First, we rewrite the generating function (23) as 



M-1 MN-1 



M-1 N-l 
2 2 



where 



Gmd(x,v) = 2Mz'"*'y""i ' JJ JJ#(m,n), 

m=l n=l 



12-2 17171 ■ 1 nn 

q(m,n)=A\z sin h sin 

M N + l 



To extend the limits of the products in (Al) to M — 1 and N as in (24), we note 



M-1 N 



|1 jj g(m,n)=C 1 C 2 ,Cz 



m=0 n=0 

(m,n)^(0,0) 



' M-1 N-l 
2 2 



-i 4 



II i i //("'•") 



m=l n=l 



M, N = odd, 



(Al) 



where C±, C2, C3 collect respective products for m = 0, n = 0, and {m 7^ 0, n = (N + l)/2}. 
Namely, for M, N odd, 

N N , s 

Ci = n^ ^) = n( 4sin2 A 7TT )=(N + l) 2 , N>1, (A2) 

n=l n=l ^ ' 

M-1 M-1 , x 

^ = n^K°)=n( 4 ^ sin2! ^)= MV(M " i) ' M>i > (^3) 

m=l m=l ^ 

(A4) 



M-1 



m=l 



Msinh -1 [ - 

z 



, M>1, 
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where the product (A4) is a special cases of the identity [19] 

Af-l 

Y[ (4 sinh 2 # + 4 sin 2 ^) = 4sinh 2 (M#), M > 1. (A5) 

m=0 

Combining these results, the generating function (Al) reduces to (25). 
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